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Abstract 



We investigate rescaling transformations for the Vlasov-Poisson and 
Euler-Poisson systems and derive in the plasma physics case Lyapunov 
functionals which can be used to analyze dispersion effects. The method 
is also used for studying the long time behaviour of the solutions and can 
J^ ' be applied to other models in kinetic theory (2-dimensional symmetric 

?— j I Vlasov-Poisson system with an external magnetic field), in fluid dynamics 

(Euler system for gases) and in quantum physics (Schrodinger-Poisson 
system, nonlinear Schrodinger equation). 
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1 Introduction 

Consider the Vlasov-Poisson system (VP) 

dtf+v-d^f-d,u-d,f = o, 



AU — ep, p{t,x)— I f{t,x,v)dv 

and the pressureless Euler-Poisson system (EP) 

dtp + div{pu) = 0, 

dtu+{u-dx)u^-dxU, 

AU^ep. 

Here i>0, x, v, u = u(t,x)€TR'^, d>l is the dimension of the physical space, 
e = +l corresponds to the stellar dynamics and e = — 1 to the plasma physics 
case. Throughout this paper, we shall assume that / is a nonnegative function 
in L°°(IR"'',L^(IR'' xIR'')). Formally, we have the following relation between 
these two systems: a pair {p,u) is a solution of (EP) if and only if 



f{t,x,v) = p{t,x)6(v — u{t,x) 



is a solution of (VP) where 6 denotes the Dirac delta distribution. In this 
situation u can be recovered from / via the identity 



p{t,x)u{t^x) — / V f{t,x,v)dv. (1-1) 

In this sense (EP) is a special case of (VP), and we will see later that the 
asymptotic behaviour of (VP) for large times is connected with a special solution 
of (EP). On a rigorous level the relation of (VP) with (EP) is investigated in 
[0. 

Throughout this paper, we assume for simplicity that the solutions of (VP) 
are of class C^ with compact support with respect to x and v, which allows us to 
perform any integration by parts without further justifications (except maybe in 
dimension 2). The results then pass to less smooth classes of solutions, assuming 
for instance that / belongs to C''(IR+,L^(IR'' x IR'')) (see for instance ^^ or 
[p9|) and is a global in time solution to the Cauchy problem corresponding to 
an initial data /o satisfying for instance: 

(d = 2) /oeL^ni°°(IEl^ xIR^) is such that for some e>0 the quantity 

fo{x,v){\x\^+' + \v\^+^ + \Uo{x)\)dxdv 



(with Uo{x) = ~^log\x\* J fodv) is bounded (see plj). 



[d = 3) /o e £^ n L°° (Si^ X H^) is such that for some e > and p>3 the quantity 

Mx,v){\v\^+' + \xndxdv 



is bounded (see |2y, g^, |2|, g^ and |g, ||] for the propagation of moments). 

For weak solutions obtained as a hmit of an approximating sequence (for in- 
stance, if we assume no moments higher than 2), the equahties have to be 
replaced by inequalities. 

For the Euler-Poisson system, we shall consider only C^ solutions. The 
results presented in this paper have to be understood as either a general method 
on how to obtain dispersion effects without taking care of the existence or the 
regularity of the solutions, or as a method to derive a priori estimates for less 
regular solutions (by passing to the limit with smooth approximating solutions). 

Our paper is organized as follows. In Section 2, we introduce linear scalings 
and explain why they give rise to singular self-similar problems. How to remedy 
this pathology with time-dependent scalings is explained in Section 3. In the 
one-dimensional case, the information on the solution is sufficient to provide 
the convergence of the rescaled solution to the asymptotic measure. Section 4 is 
concerned with the Lyapunov functionals and constitutes the heart of this paper: 
the energy of the rescaled system turns out to be a Lyapunov functional for the 
initial problem. A more straightforward (than the full time-dependent scaling 
method) approach to the Lyapunov functionals is also given. In Section 5, we 
use the Lyapunov functionals to describe the asymptotic behaviour (dispersion 
rate) of the solutions in the plasma physics case. 

Rescalings for the study of large time behaviour have been widely used in 
various fields of applied mathematics but appear to be rather new in the context 
of kinetic equations: in that direction we may mention the studies made by 
J. R. Burgan, M. R. Feix, E. Fijalkow, and A. Munier (see [Q) and J. Batt, 
M. Kunze, G. Rein in [nl. Our main point is to make the link between rescalings 
preserving the L^-norm and Lyapunov functionals (or pseudo-conformal laws) 
and to explain on various examples of conservative systems why it actually 
provides a general method for the study of large time asymptotics. 

While Sections 1-5 are exclusively devoted to the Vlasov-Poisson and Euler- 
Poisson systems. Sections 6-8 are concerned with other problems of kinetic 
theory, fluid mechanics and quantum physics. The relation between these var- 
ious domains has been noticed for a long time (see for instance [p6|), but it is 
surprising that the estimates given in |1^ have been adapted to kinetic mod- 
els only recently. Here we proceed in the reverse historical order, from kinetic 
equations to fluids and quantum physics, and this approach actually seems to 
be very powerful. 

To conclude with the introduction, it is worth mentioning that many of the 
estimates we are giving in this paper were already at least partially known. The 



point is that we present a systematic and elementary method which takes the 
nonhnearity of the model very well into account (this was not necessarily the case 
in the preceding papers) and gives rise to a more precise form of the Lyapunov 
functionals (in the sense that these Lyapunov functionals also include second- 
moments in the x-variable) which are natural for the problems we consider. 

2 Linear scalings 

Let f = f{t,x,v) be a solution of (VP). Then for any X,fi>0 

is again a solution of (VP), 

3t/A,/i + V ■ dxfx,t, - 9:eC/a,m • dvf\f, = 0, 



AC/a,p = £Pa,m, PA,M(i,a:)= / f\^f,{t,x,v) 

JM'' 



. ,.,^ , , ,dv 



with 



P\,p.{t,x) = X p{Xt,nx), 

Uxjt,x) = x^^i-^u{xt,^JLx), 
d^UxJt,x) = x^^i-^^.,u{xt,^lx), 

as can be checked by direct computation using for instance the following integral 
representation of dxU: 



dxU{t,x)^- * I fit,x,v)dv. 

Here S"'^^CIR is the unit sphere of dimension d—1. Similarly, if {p,u) = 
{p,u){t,x) is a solution of (EP) then for any A,/i>0 

P\,p.{t,x) = X^p{Xt,px), ux^^{t,x) = Xp^'^u{Xt,px) 

is again a solution of (EP), and the potential is transformed as for (VP). This 



also follows from Relation (1.1) between (VP) and (EP). 

If we require that the L^-norm of p(i), which is a conserved quantity for 
(EP) as well as for (VP), is preserved by the scaling, A and p must satisfy 

X^p-'' = 1 

and the rescaled distribution function is fx x^/-^- -^ standard way of studying 
the asymptotic behaviour of / would then be to consider a self-similar solution, 
i.e. a solution which satisfies 

f{t,x,v) = h^.,,{t,x,v)^X^-''f{Xt,X^/^x,X^/''-^v) 



for any A>0. This solution would then be given by its self-similar profile 
f{^,ri)~f{l,£^,ri) (choose A to be j). Then 

is a solution of (VP) if, at least formally, / is a solution of 

in the new variables ^ — t^dx, rj — t^^dy. However, it is clear that as t^0+, 
g{t,x,v) does not converge to a well defined measure for which one might estab- 
lish an existence result, except for d = l. This difficulty is completely removed 
by considering general, non-singular, time-dependent scalings. 

3 Time-dependent scalings 

Consider the following transformation of variables in (VP), where the positive 
functions A(t), R{t), G{t) will be determined later: 

dt = A^{t)dT, x^R{t)£,. 

Thus, assuming that t^^x{t) and t^^S,{t) satisfy ^ — v and -^—f] respectively, 
the new velocity variable 77 has to satisfy 

Let F be the rescaled distribution function: 

f{t,x,v)^G{t)F{T,C,v)- 

The aim is to choose this transformation in such a way that the transformed 
Vlasov equation is still a transport equation on phase space and contains a 
given, external force and a friction term. The inverse transformation is 

dT = A-\t)dt, ^ = R-\t)x, ry=.^'^^^ ^- ^^^^^ 



R{t) V R{t) 

Here ' always denotes derivative with respect to t. If i' and W are defined as 
the rescaled spatial density and the rescaled potential respectively, then 

v{r,0^j /ir,(,v)dv^j^p{t,x), 



A 



2d 



A 



2d 



and the Vlasov equation transforms into 

- R^£,-d,^F-R'^GA'^-'^'^d^W-drjF + A^%F = Q. 
R G 

We want this to be a conservation law on (^,7/)-space, so we require 

A i?_ 1 G 
A~R~2dG 



which holds if and only if 



-"'^ 



2d 



(3.1) 



(3.2) 



for some constant c > 0; recall that G should be positive. For simplicity and 
without any loss of generality, we may choose c=l, and the Vlasov equation 
becomes 



drF + r]-d^F + diY^ 



a G R I 



Next we require that the external force in the above Vlasov equation becomes 
time-independent and that there is no time-dependent factor in front of the 
nonlinear term, i.e. 

it—— = — £Co, 



R 



(3.3) 
(3.4) 



where cq > is an arbitrary constant. In view of (3.2) we get 

and R has to solve 

R + ecoR^-'^^O. 



(3.5) 



At this point, we may state the following remark: 

Remark 3.1 Every solution of Equation ^3.^ ) has the following properties: 
(i) For any A>0, ti-^ Rx{t)^CQ '^X^dR[Xt) is a solution of 

R + eR^-'^ = 0. (3.6) 

Without loss of generality we therefore assume that cq = 1 in what follows. 
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(ii) With Ro = R{0) and Rq = i?(0) we get, for d^l, 



R{t)- 



-t +RQt + RQ 



If d>2, it is easy to carry out one integration of Equation (S.t): 



^R\t)+e\ogR{t) = ^Rl+e\ogRo for d^2, 



5«'W-5^ 



R 



2-d 



1 



it) = oRo--r^Ro ford>3 



2"u ^_2 " 

In the plasma physics case, R(t) cannot change sign and is well defined 
for any t G TR. Moreover, 



\ogR{t)^-R\t) 



ii?2 + logi?o>-^i?o + logi?o ford^2, 



Q-CR^^'^it) 



d-2 



Rn 



ie 



i?, 



2-rf 



< 



d-2 



-R(\ 



-Rl^"^ ford>3. 



Together with Equation (3.6) this proves that there exists a unique tQ(£TR 
such that -R(to) >0 and R{tQ) — 0, and R{t) > Rq for any t^to, provided 
e = -l. 



(Hi) If tis- R{t) is a solution of Equation (S.t) with e = —l, ti-^R{t-\-a) is 
a solution too for any given a G ]R. Combining this with the invariance 
through the rescaling Xi-^ R\{t) with R\{t) — X^d R(^Xt), we may always 
require Ro = l and Rq — O without loss of generality as long as we are 
interested in the asymptotic behaviour of f when t—^+oo. Note that with 
this special choice for Rq and Rq, at t — 0, G(0) = j4(0) = 1, and if we 
assume r(0) — 0, then 

S,(t = 0,x) =x, ri{T = 0,x,v) —V and f{t — 0,x,v) — F{T — 0,x,v). 

The time- dependent rescaling has the interesting property that it does not 
introduce any singularity at t — 0, and with Rq — 1 and Rq = 0, the initial 
data for f and F are the same. 



(iv) The singular self-similar solution (2.1) corresponding to the linear scalings 
of Section 2 is — when it exists — the solution one expects to get in the limit 
case Ro = 0. Formally, this solution also corresponds to the limit of R\[t) 
as X^ +00 . 



(v) For e = —1 and d>l, -S^t as t^ +oo and 

R{t)^t^ 

R{t)r^t 



for d—1, 
ford^2, 
ford>3. 



(3.7) 



With R solving Equation ( |3.6| ), we obtain the foUowing rescaled Vlasov- 
Poisson system (RVP): 



drF + r]-d^F + div,^ 



e(-d^W+^Ri-^Rrj]F 



= 0, 






The relation between the old and the new variables is 

dt = R'^/^dT, dT = R-'^^^dt, 
x==:R^, ^^R-^x, 
v = R^ + R^ ^1], ri = R^ Mw-— X 
and the rescaled functions are given by 

W{T,()^R'^-^U{t,x), d(^W{T,i)^R'^-^d^U{t,x). 
If we consider (EP) we find that with 

v{T,O^R'^^^[u{t,x)^-x 

the rescaled Euler-Poisson system (REP) is 

dTV + d\v{vri) = 0, 

d-4 



drri+{r]-d^)r] — e^ — d^W' 



-R^-^Rf], 



Note that this rescaling as the one for the Vlasov-Poisson system introduces a 
harmonic force term e^ and a friction term which is proportional to the velocity. 

There exists a unique steady state with a given L^-norm M for which the 
particles are at rest and uniformly distributed in the unit ball centered at 0, 
and the self-consistent force is exactly balanced by the external force. Define 



F^i^,v) = ^ZiOKv) 



where S is the usual Dirac distribution. If F^ and (^^,^,'7;^ = 0) are the sta- 
tionary solutions of (RVP) and (REP) respectively such that 

ll^oo \\l^{TR''xTR'') = lkoollLi(]R") =M, 

then 
and 

^oo(C)=C?- lBd((ji//|5d-i|)i/d). (3.9) 

Here B'^{r) denotes the ball with radius r centered at SIR , and I^ denotes 
the characteristic function of the set w. The inverse rescaling transformation 
takes this steady state into 

f^{t,x,v) = -^j7Y'^B-^{R{i)(M/\s-^-^)^/-^){x)5\^--—x\ , (3.10) 



and 



A'f / \ fT / \ A4 f \ '^\ ) 



It is easy to see that this defines a weak solution of (VP) or (EP) respectively. 

In the plasma physics case we have R{t) > for t > provided i?(0) > so 
that for d < 3 the particles are slowed down by a friction force, and on physical 
grounds one would expect that the steady state written above is a global at- 
tractor for (RVP) or (REP) respectively. In [Q this was carried out rigorously 
for the case d=l. We will see in the next section that this is not true in gen- 
eral, at least in dimension d = 3 for (EP). However, the rescaling still provides 
informations on the asymptotic behaviour of the original system for large times: 
the energy for the rescaled system gives rise to a Lyapunov functional for the 
original system by which dispersion effects and the asymptotic behaviour can 
be analyzed. 

4 Lyapunov functionals 

In this section, we investigate the behaviour of the total energy of (RVP) and 
(REP). Let us consider first the case d>3. The potential energy term is the 
same for both systems, namely 

For (RVP) the kinetic energy reads 

EUr)= ff\v\'Fir,^,v)dvdC, 



while for (REP) it reads 



Sk(r) = y'|77p(r,0^^(r,Ode 



Recalling the remark on the relation between (VP) and (EP) from the intro- 
duction, the second formula can be viewed as a special case of the first one, and 
for both systems we find after a standard computation: 



drV 



(E^{T) + Ep{T)]={d~4)Ri-'^REi,{T) for d>3. 



(4.1) 



This result is easily achieved as follows. First integrate the rescaled Vlasov 
equation with respect to 77, 



i9r'^ + div(;/77) = 0, 



(4.2) 



which is nothing else than the local conservation of mass. Then multiply the 
rescaled Vlasov equation by (|?7p— e|CP) ^^nd integrate with respect to ^ and 77, 
and then by parts. 



i\v\'-e\^\')ivdiF + e^-d,F)d^dv^O, 



so 



d7//^l^|- 



e\S.\^)Fd^dTj + 2 // 7j-d^WFd^drj^{d-A)R 



^-'RJJH 



■Fd^dr]. 



An integration by parts with respect to f gives 

ffrj-d^WFd^dr]^- fw{d^- fr]FdT])d^, 



which using (4^) gives (4J). Note that (12) is written for instance in the sense 
of distributions and that in the above integration by parts one has to check that 
no boundary term appears. This is true for d>3, but not for d~2 aswe shall 
see below. 

RecaU also that R'^/^R^ R'^^^dR/dt^d[R{t{T))]/dT. Let us rewrite the en- 
ergy for the rescaled systems in terms of the original variables: if we define P 
and K by 



P{t)^E^{T{t)) = R'^-\t) \U{t,x)~e 



RHt) 



p{t,x)dx, 



and 



K{t) = E^{T{t)) = R''-\t) 



R 

V X 

R 



f{t,x,v)dxdv 
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for (VP) or, for (EP) 



K(t) = E^iT{t)) = R^-'{t) 



u{t,x) X 



p{t,x)dx , 



then because of (4.1) 



L{t)^K{t) + P{t) (4.3) 

is a non-increasing quantity with respect to t for d — 3, 4: 

§ = «,-4,|«<0. (4.4) 

Because of the integrations by parts in the intermediate computations, the above 
formulas are true only for d>3. We will now consider the cases d=l and d = 2. 
In dimension d~l with e — —l (plasma physics case), direct computations 
involving the kinetic energy and integral quantities related to the force field 
have been used in pi to prove the exponential convergence (in the rescaled time 
variable r) of F(t,-,-) towards Foo in {W^'°°(iR^)y and of d^W{T,-) towards 
ds^Woo in L^(Si). The same computation also holds true for the solution of 
(EP) if it exists globally in time: 

Proposition 4.1 Assume that d—1, e — ^1 and consider a global solution 
(t,x)^^ {p{t,x),u{t,x)) of (EP) in C^(IR+x]R) such that for any t>0, p{t,-) 
has a compact support. Then 

Hr{t),o - R{t)p{t,R{t)o , v{T{t),o - -i^ (^u{t,R{m ~ Rm 

with T(t) = 21og(l + i) and R{t) = {l + t)'^ is a s olution of (REP) and converges 
to {v^,0) where v^ is given by Equation 4?. 4 ), with Af = ||p(i,-)||ii(]f;) ; there 
exists a positive constant C such that 

while the electric field ~d(W{T^-) — J_^u[TX)dC, — j\\v{T,-)\\i^ifjg_) converges in 



L^(IR) to ^d(W^ which is given by Equation (3.8): 

In terms of the original variables and with the notation of Section 3, this means: 



[i+ty p{t,{i+tY-)-v'^\\i^w^.^-^ 



< 



c 



\d^u{t,{i+tY-)-d^w^\\L^^)< 



_c__ 
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Note that in Proposition 4.1, we made for R{t) the same choice as in |l|, 



which means that with the notation of Remark |3.1| we consider the solution of 
Equation ( |3.5[ ) corresponding to i?o = 1 and Ro = 2. 

The proof essentially follows the same arguments as in M. 

In the case e = +l (gravitational case), essentially nothing is known con- 
cerning the asymptotic behavior of the solution. If d— 2,3,4 and e = —l, the 
question of identifying the limit of F(t, •, •) or j^(t, •) in the sense of measures as 
t^Too = Jq °° R~'^^'^{t)dt (which is finite as soon as d>3) is an open question. 
As already noted, a natural conjecture would be to identify this limit with Fjjf 
for the solution of (RVP) and i/^ for the solution of (REP) as in dimension 
d=l. In other terms, the stationary state of the rescaled equation would be an 
attractor for the solutions of the rescaled system in dimension d>l. li d>3, 
this is not true in general. 

Counter-examples. Consider a solution for which it is the case and shift 
the initial data by a constant velocity. Since asymptotically the support of the 
unsealed solution grows linearly in time, after rescaling, the shifted solution 
cannot converge to the stationary profile. One may then ask the same question 
in the reference frame of the center of mass. The following counter-example for 
(EP) again shows that for d>3, e = —l, the answer is negative. 

Consider in IR the solution corresponding to the following initial data: 

p{t = Q,x) = 3Ib3(i) [x) + Ib3(3)\b3(2) (x) , 
u(i = 0,a;) = if \x\<l, 
u{t = 0,x)^x if 2<|a;|<3. 

For any t>0, the solution is supported in the union of a centered ball of radius 
R{t) (which obeys to Equation ( |3.5D J and of a centered annulus of inner radius 
Ri{t). A straightforward computation shows that R and Ri satisfy 



R^-^, R{0) = 1, RiO) = 0, i?i = ^ , i?i (0) = 2, i?i (0) = 2 , 



^,i?(0) = l,7?(0) = 0, i?i = -^2 
respectively, and an integration with respect to t gives 



i?2(i) = 2-^-<2<4<5 %-^R^(t) 

^ ' R{t) Ri{t) 1^ ' 

for any t>0. As t^ +00, y/2 = lim^+00 —j^ < lim^+oo -^^ = V^ which again 
forbids the convergence to the stationary solution after rescaling. 

In dimension d = 2 for e = — 1, the situation is different: ti—fR{t) grows su- 



perlinearly (as for d=l: see Remark 3.1, (v)), and the question is still open 
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Consider now the case d = 2 for (VP). The main difSculty comes from the 
integration by parts, and one has to be very careful with the terms involving the 
self-consistent potential U since VU essentially decays like l/\x\. Let {p,pu) = 
Jfit,x,v){l,v)dv, {iy,iyr]) = J F{T,^,r]){l,r])d'i] and M= ||/(i,-,-)||ii(]R2x]R2). 



2-H){r,0d^ 



Rx 



2 R^p{t,x){u{t,x) — —x) 



R dx 



-R I p{t,x)dx + R 
1 



• (pu){t,x)dx 



'MR~R 



dtp{t,x) dx 



using the local conservation of mass dtp + dx{pu) = 0. Similarly. 
'd^W-{friF{T,^rj)dr,)d^ 

M^ R 1 ^ / 

pU{t,x) dx 



— — - / U{t,x)dtpit,x)dx 



AP R 1 d 
^ '2nR^2dt^ 

271 R 2dT J '^ '^' ^' 

cf. [Q for more details. Of course the computations are exactly the same for 
(EP). Thus in dimension d=2, the definition (|4.3|) has to be replaced by 



L{t)^K{t)+P{t) + — \ogR{t), 



(4.5) 



so that Equation (4.4) still holds, see also Remark 4.4 



Equation (4.4) provides an identity which is a sharpened form of the Lya- 
punov functional (also called pseudo-conformal law: see Section 10 for the re- 
lation with the Schrodinger equation). A simple form of this identity had been 
discovered independently by R. lUner and G. Rein, and by B. Perthame, cf. 
[ |l6| , |25|| . The improved Lyapunov functional has the striking property that it 
easily provides all the terms that one has to take into account in the case d=2 
(see |ll| for (VP)) in a quite straightforward manner, while a direct approach 
was far from being obvious. 

Theorem 4.2 Assume that f is a s oluti on of (VP) with AI ^\\f{t,-,-)\\ and that 
if— >i?(t) is the solution of Equation (S.t) with R[Q)~\, i?(0) = 0. The function 
ii— > L{t) given by 



L{t) = R 



d-2 



M^ X R"* 



R 



fdvdx + R 



d-2 



W' 



U{t,x)-e-^ ) Pdx 
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for d>3 and 

Lit): 



]R2x]R2 



R 

V X 

R 



f dv dx - 



E,^ 



U{t,x)-e^] pdx+—\ogR 



if d — 2 is decreasing for d~2,3, constant for d — A, and for any d>2 satisfies 



dt ^ ^ 



R 

V X 

R 



fdvdx. 



Moreover in the plasma physics case e~—l, L is bounded from below, and for 
d^2,'S, 

2 



R{s)R''-'{s) 



R^'xIR'' 



R(s 

V -^x 

Ris) 



f(s,x,v)dvdx j ds < +00 



Proof. dL/dt has already been computed above. For e = — 1, the proof of the ex- 
istence of a lower bound is straightforward except maybe for d = 2. In that case, 
the Lyapunov functional is decreasing but might a priori be unbounded from 
below, and we have to estimate it. This can be done using M = J-^2p{t,x)dx 
and Jj^2pit,x)\x\^ dx by splitting the integral with respect to x and y into two 
parts corresponding to \x — y\<k and \x — y\>k, but a more straightforward 
approach can be deduced from Jensen's inequality using the fact that —log is a 
convex function: 



2ttJJj, 



H-'xIR^ 



\og\x - y\ p{t,x)p{t,y) dxdy 

m2 rr 

47r iiiR^xii- 

47r 
An 



> 



> 



dxdy 

M2 

i2 / \ / \ dxdy\ 
x-y\ ■p{t,x)p{t,y)-—^ 

]R2x]R2 M ) 



-log|a;-yp ) ■ p{t,x)p{t,y) 
■x]r2 \ 

log 

log(2//M) , 



where I = J\x\^pdx, and an optimization on / > gives 
M2, „ I 

which proves the result. □ 



— log(2//M) > ^ [1 - log(M/2^)] , 



Remark 4.3 Fors = +l andd = 3 or A, using the Hardy-Littlewood-Sobolev in- 
equality and classical interpolation identities, one proves that the self- consistent 
potential energy term |lVt/||^2 is bounded in terms of K by 



IVC/II 



<rr'\\ f||2(l-(<i^-4)/4d)|| ru(d-2)/d j^ (,x(d-2)/2 

^d) !i<^||/||rlfn^d„n^d-, I U 1 1 Too /TRd v TRd-l ^^ V^-* ' 



\L^{TR.^) ^^lU llLi(]R<ixIR'^) II-' " L°° {TRf^ xW^ 

see Section 5 for more details on interpolations. 



14 



Remark 4.4 In dimension d — 2, for e~—\, it is probably easier to compute 
dL/dt and prove (4-4) directly from (VP) using the identity 

f f [ NP 

If {x-v){dxU ■dvf)dvdx = — {x-dxU)pdx—— — 

once the equation for R is known, cf. |7j|/. 



Note that with the help of (3.7) and the resuhs of Theorem 4.2, we recover 
the results of IK , 25 in dimension d = 3 as well as the results of [[ll[ in dimension 
d = 2. Very similar results of course hold for (EP) since the estimates on the 
Lyapunov functional in dimension d = 2,3,4 are the same. 

Theorem 4.5 Assume that {p,u) is a global strong solution of ( EP) with M = 
||p(t, •)||^i(][^d) and that ti~^R{t) is the solution of Equation ^3.f\) with i?(0) = 1 
andR{0) — 0. The function ti-^ L{t) given by 



Lit) = R'^-^ f 
for d>3, and 

Lit) 



R 



pdx + R 



d-2 



R" 



U-e- 



Rd 



pdx 



]R2 



R 



\xV\ ]VP 

pdx+ I lU-e'-^] pdx+ --logR 
f]R2 V R"^ J 2n 



for d — 2, is decreasing for d — 2,3, constant for d = A, and for any d>2 satisfies 

2 



'§{t)^id~A)RR'-' 

ox ./ ]R ° 



u{t,x) X 

^' ' R 



p{t,x)dx . 



Moreover in the plasma physics case e — —l, L is bounded from below, and for 
d = 2,3 



+ 00 

R{s)R'^-^{s) 

yJTR" 



R{s) 
R{s) 



p{s,x)dx I ds < +CX) . 



The case d — A appears to be the limit case to which the above method for 
finding Lyapunov functionals in the plasma physics case applies since for (i> 5, 
ti—>L{t) is increasing. However for d>4, we may write 



2 L 



<0, 
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and thus obtain 



TR'' X IR,'' 



R 



f{t,x,v)dvdx = O (R-^) = O (f-2) 



since for d>3 all the quantities involved in L{t) are nonnegative and R{t)^t 
as i^ +00. 

In this last part of Section 4, we will derive the Lyapunov functionals in 
another way, not because of the case d>4 (which is of minor interest for (EP) 
or (VP) in itself), but because the method is simpler and will be applied to other 
systems in Sections 6-8. We assume that e — —l in the rest of this section. 

We may indeed notice that all the quantities we have been taking into ac- 
count are integrated in the x variable, so that the change of variable ^(i,a;) = 
x/R{t) does not play any role in the estimates. Let us first consider the Vlasov- 
Poisson system (VP). According to the above remark, we may use the change 
of variables 

ri{t,x,v)=v-—x, f{t,x,v)=:^F{t,x,T]) 
so that F solves the rescaled system (R'VP): 

dtF + 7j-d,F-^x-dr,F-d,U{t,x)-dr,F+^[d^{xF)-dr,{rjF)]^0, 



'd,U{t,x) 



\S''-^\\x\'i J J,. 
As for (RVP), we may compute the energy 



F{t,x,ri) drj . 



E(t) = 



H^'xIR'' 



-I 



U]Fdxdri 



if d > 3 and 



E{t)^ 



r^\^ + ^\,\^ + u\Fdxdr^+^^ 



' R'' X R' 

for d = 2. This energy is a decaying function of t : for any d>2 



dE 
Itt 



{t) = -{d-2) 



R'' X E,'' 



FU dxdrj 



d , R , ^R R 

— ( — ) + 2 lla; 

/ TO d y TR 1^ \ -^ -^ -^ 

We may now define L{t) ~B{t)E{t). For any d> 3, 



^-2|m^' 



F dxdrj . 



dL 



B-{d-2)^B 



wurdx 



n-^ 
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B-2^b) ff F{t,x,r^M^ dxdij (4.6) 

^l+f^(|) + 4l)^) // F{t,x,v)\x\'dxd^, 

n \at K UK/ / J Jm''xM'' 

while for d—2, 

— = (mogR-B-)M^ 

+ (b-2^b)[[ F{t,x,rjM^dxdv (4.7) 

H \at K UHJ /JJ]R2xIR2 

For d>3, the following conditions are sufficient for L to be nonincreasing: 

1) B{t)^R{tY-^, which implies B- {d-2)BR/R<0, 

2) d<4, which implies B~ BR/R^-{4:-d)BR/R<0 , 

3) i? = i?P, i?(0) = l, i?(0) = with p<-((i-l), which implies (^(|) + 



and we recover the results of Theorem (4.2) for (i = 3,4; for d = 2 take B — 1 



Remark 4.6 If d>2 (including the case d>A), we may choose B = R'^~^~^, 
0> max(0,d— 4), and R solving the equation R = Rp, i?(0) = l, i?(0) = for 
some p<9 ~{d~l) without any further restriction on d. Note that for (i>4, 
p<—l and <d—2, one recovers the estimate one would have for the free trans- 
port dtf +v-dxf =0, since in that case f{t,x,v) = fo(x — vt,v) and 

f{t,x,v)\x-'Vt\'^dxdv— // fo{x,v)\x\'^ dxdv. 

For the consequences on the dispersion rate, see Section 5. 

An analogous method also works for the Euler-Poisson system (EP). If we 
shift the velocity u{t,x) by an unknown "bulk" velocity ^x, so that ri{t,x) = 
u{t,x) — ^x, then {p{t,x),ri{t,x)) solves the system (R'EP): 

dtTl+j^{-)x+{{Ti+-x)-d^)r^+-{7^+-x) = ~d^U{t,x), 
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As for (RVP), we may consider the energy: 

E{t)^J(\r,it,x)\'' + ^\x\^ + U{t,xy\p{t,x)dx 

for d>3 (if d~2, one has to add the term %f^)- Again tt-^ E{t) is decaying: 
for any d>2, 



f'-*"-^*!!/^"- 



w 



d R R R 
di^R'^ RR 



4h 



pdx. 



As for (VP), we may also define L{t) = B{t)E{t), and the rest of the discussion 
is exactly the same. 

The method for finding a Lyapunov functional can be summarized as fol- 
lows: first change the velocity variable by subtracting a velocity -Sx for some 
increasing function R, then compute the energy associated to the new equation 
and finally choose the Lyapunov functional to be L{t) — B{t)E{t) where Bit) is 
the function of t which has the maximal growth in order that L{t) is still a de- 
caying function of t and corresponds to a function ti— > R{t) solving an adequate 
ordinary differential equation which takes the nonlinearity into account and has 
to be chosen well. Of course, one way to find an equation for R is to apply 
the method of the time-dependent rescalings of the beginning of this section. 
This method is sufficient to extract the asymptotic rate of decay of the relevant 
quantities, as we shall see later in several other cases, cf. Sections 6-8. 

5 Asymptotic behaviour, dispersion 

An estimate of the rate of dispersion of a solution / of the Vlasov-Poisson system 
(VP) in the plasma physics case e — —l is given by the interpolation of p{t,x) = 
J-^a f{t,x,v) dv between the _L°°-norm of /, which is preserved for strong solu- 
tions, and the momentum JJ^a^^j^af\v—j\^dxdv^JJ^dy^^df\v~^x\'^dxdv as 
t -^ +CX3: there exists a constant C = C{d) > such that 



JR"* 



fdv 



<C- 



Ld + 2(]Rd) 






TR'^xTR.'' 



R 



dxdv 



(5.1) 
for a systematic study of these interpolation inequalities see [|ll| and references 
therein. 

The asymptotic form of the Lyapunov functional was given in p6| , p5[ for 
the case d — 3 and in |1^ for the case d — 2. Using R{t), we remove the difficulty 
due to the singularity at t = and recover the known results. The use of the 
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decay term of the Lyapunov functional allows us to prove that the decay is not 
optimal. 

Proposition 5.1 Assume that f is a strong solution of (VP) in the plasma 
physics case £ = —1 and ti-^R{t) is the solution of R = R^^'^ with R{Q) — \ and 
-R(0) = 0. Then f obeys to the following Strichartz type estimate: if d^2 or 3, 



+ OC 



R''-\t)R{t)[ // f{t,x,v) 



R{t) 

V -—X 

R{t) 



dxdv\dt<C, (5.2) 



and for d — 3,4:, we have the following dispersion estimate 



\\p{t,-)L^,^,<CR{t)-''—^ 

L " [ax ) 



, a — z 



(5.3) 



Here C denotes various positive constants which depend only on d and fo, and 
L is the Lyapunov functional of Theorem {A. Ifd — 2, 



liminf ||p(t)||i2(]R2)=0 



t^-\-oo 



(5.4) 



The proof follows from Theorem [4.2| and the interpolation identity (5.1) 
given above. The decay of p{t,-) in L^(JR^) is given by the decay term of the 
Lyapunov functional; see Remark |5.2| below. Estimate ( |5.2| ) for d = 2 has been 
improved compared to [ [ll| . 

Remark 5.2 The decay given in Proposition |5.j| is not optimal. Consider in- 
deed a function ti—th{t) such that h>\, ]iTat^j^ooh{t)~+oo and 



+ OC 



■ — +CX3 if d — 2 , 



+ OC 



ds 



his) 



~ +00 if d = 3 . 



JO sh{s) 
For instance, one may take for t > 

h{t)^log{t + 2) if d=2 and h{t)=tlog(t + 2) ifd='d. 



The hounds {5.i) immediately provide for d~2,2i 



liminfft,(i) // f{t,x,v) 

t^+oo JJu-'xTR-' 



R{t) 

V -—X 

Rit) 



dxdv = 0, 



and as a consequence, the decays in (5.5) and (5.4) are not optimal 



Similar results can of course be obtained for any d>5, using Remark 4.6 



We may notice that the decay in (5.3) for d> 4 is the one which is obtained for 
the free transport equation when considering the second moment in a; — vt. 
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Remark 5.3 The Lyapunov functionals given in / pj| , \Tq, |^q / correspond to the 
asymptotic form of R(t) as t— ++00. The fact that this asymptotic form also 
gives a Lyapunov functional is easily explained by the scaling invariance of the 



equation (see Remark 3.1): if one replaces R{t) by t for (i = 3,4 or t^/\ogt for 



d—2 in the expression of the Lyapunov functional L{t) of Theorem l^.i., L{t) 
would still be a Lyapunov functional (consider its limit as X—>-+oo). 



Similar results for the pressureless Euler-Poisson system (EP) also hold ex- 
cept that no direct interpolation can be used. The decay only holds in a weak 
norm defined as follows (assume here that d>3): let us consider the space 

pi'2(]R'^) = |0eL#2(lRrf) : V0eL2(]R'*)|. 
and define on its dual space the norm 

IIIpIII = ||p||(pi,2(R.)), =sup|y p0da: icj^ev'^^iiR"), |1V0|L.(R.) < i| . 

If U&V^''^{]R'') is such that -AU = p, then |||p||| < ||VC/||i2(]R,d). Using the 
same notation as in Section 4, if (i = 3,4, there exists a positive constant C such 
that 



p{t,x)\r]{t,x)\'' dx^ R{tf-^ / p{t,x) 



. . R 

u(t,x) X 

^ ' R 



dx<C 



and 



Rit) 



d-2 



pit,x)U{t,x)dx:^R{tf-^ / \VU{t,x)\^dx<C. 

TR.'' JtR'^ 



The last inequality can be reinterpreted as an estimate on the weak norm 1 1 1 - 1 1 1 
oipit,-). 

Proposition 5.4 Assume that {p,u) is a C^ solution on M"'' x M of (EP) in 
the plasma physics case e = —l and ti— >i?(i) is the solution of R + eR^~'^ = 
with R{0) — 1 and R{0) — 0. Then {p,u) obeys to the following Strichartz type 
estimate for d — 2, 3; 



+00 



R''-'\t)R{t) 



Moreover, if d — 3,4, then 



p{t,x)\u{t,x) — 'j^TT^l dx ] dt < +00 . 



limsupi2-i|||p(i^.)|||=limsupi?(i)2 ^\\\/U{t,-)\\^2(^^d) < +00 . 
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Remark 5.5 Ifd~2, we cannot use the {V^'^iJR ))' -norm as in the case d>3, 
but the following estimates for the solutions in the plasma physics case e — —l 
of the pressureless Euler-Poisson system (EP) hold: 

1 /■ M^ 

lim -, T^TT I p(t,x)U(t,x)dx = —— — , 

t->+cx) [ogK[t) Jjn2 2tt 



1 /■ / M ,2 , M^ 



/R2 

These estimates are easily deduced from the conservation of the energy, the 
expression of the Lyapunov functional L(t) and the estimate given in the proof 



of Theorem ^.i 



Maybe more interesting is the observation (see [|ll|) that for d = 2, which 
is the hmit case for dispersion results, the dispersion estimate gives a lower 
bound for the growth of the support of a solution corresponding to a compactly 
supported initial datum: 

Corollary 5.6 Consider for d = 2 solutions of (VP) or (EP) corresponding to 
compactly supported initial data. Assume that r{t) is the minimal radius of the 
balls containing the support of p(t,-). Then there exists a constant C>0 such 
that 

r{t)>CR{t) as t^+oo. 

Proof. As in [|ll| one may simply notice that 

^ (logRit) - log(2r(i)) j < L{t) < i(0) . D 

6 The 2-diniensional symmetric Vlasov-Poisson 
system with an external magnetic field 

In dimension d—2, we may consider the following system (VPM) 
dtf + v-d,f+ (~d:,Uit,x) + BovA ■dyf = 

X 

corresponding to a system of particles with a self-interaction through electro- 
static forces, in the presence of an external constant magnetic field Bq. Here we 
use the notation 

V2) \ Vi 
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-d^U{t,x)= * / f{t,x,v)dv 



For the linear system without self-consistent electrostatic forces, all the charac- 
teristics are circles and a solution with an initially compact support will remain 
supported in a fixed compact set for all time. With a self-consistent Poisson 
term, the situation is radically different since we get the same estimates as for 
the Vlasov-Poisson system without a magnetic field. 

We may indeed shift the velocity variable ri{t,x,v) —v— ^x, and the new 
distribution function f{t,x,v) = F{t^x,rj) obeys to the system 



R 



+ Jl[dAxF)-d,,{^F)] ==0, 



-d^U{t,x) 



27r|x|2 



F{t,x,ri) drj . 



R^ 



Exactly as in Section 4, the energy is 



E{t)^^\ogR+ ff F{t,x,r^)[\rj\' +^-^\x\'' + U{t,x) ] drjdx 



|2 I -"-1^12 



and decays according to 



dE 
'dt 



R^'xIR'' 



F{t,x,r]) 



^R, ,9 fd.R. ^RR\, ,9 „R, 1 



dxdrj . 



If / is radially symmetric, i.e. depends only on i, \x\, {x ■ v) and \x\'^ \v\'^ — {x ■ v)"^ , 
then the analogous property holds for F: F only depends on i, |a;|, (x-rj) and 
|xp|?7p — (x-?7)^, and 



x^ ■ I riF{t,x,ri) dr] ] dx= // (x-ry )F{t,x,r]) — . 

]R^ \ Jr.^ ) J Jm'^xM'^ 

The system has the same Lyapunov functional as (VP), and we obtain the same 
dispersion results as for the Vlasov-Poisson system: 

Proposition 6.1 Let d = 2. Assume that f is a solution of (VPM) and that 
ti-^R{t) is the solution ofR—-^ with R{Q) = 1, _R(0) = 0. The function ti-^ L{t) 
given by 



L{t) = ^-\ogR- 

ZTT 



R 

V X 

R 



is decreasing, bounded from below and satisfies 

R 



f{t,x,v)dvdx+ iu+-^] pdx 



dL 
'dt 



R 



H^xH^ 



R 

V X 

R 



f{t,x,v) dvdx . 
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Moreover 



+00 






]R2x]R2 



R{s) 
R{s) 



f{s,x,v)dvdx j ds < +00 



liminf ||p(t)||z,2(]R2)=0. 



7 The isentropic Euler system for perfect gases 

As another example, which does not belong to the field of kinetic equations, we 
consider the isentropic Euler system (IE) for perfect gases (for 7> 1) 

9tp + 9a;(/9M)==0, 

dtu +{u-dx)u = -dxP , 



The method goes exactly as for the pressureless Euler-Poisson system (here 
we use the second method of Section 4): the rescaled system (R'lE) given by 
ri{t,x) — u{t,x) — ^x is 

dtp + d^{p{'i]+—x)) = 0, 

The last equation can be rewritten as 

dtri + ri-dx'q+—x-dxri+—x+—ri = ~dxP^^^ ■ 
H, li Jri 

If we define the energy by 

£;(t)= / p{t,x)\r]{t,x)\'^dx+— / p{t,x)\x\'^dx+- / p^{t,x)dx, 

a Lyapunov functional is easily exhibited by considering L(t) = B{t)E{t). The 
energy is indeed decreasing: 



dE 
'dt 



M" 



„R, ,9 fd.R. „RR 



7 R 



dx 
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so that 



dL 






p{t,x)\r]{t,x)\'^ dx 



TR'' 

d , R . R R , ^ 
dt^R^ + ^RR^'' 



p{t,x)\x\'^ dx 



--fs-(7-l)di?| 



p'^{t,x)dx, 



and sufRcient conditions for L to be decreasing are therefore given by: 



I) B^Ri with g<min(2,(7- 
l)BR/R<Q. 



l)(i), which impUes B 



2|B<0 and^-d(7- 



2) R:^RP with p< -((7+1), which imphcs Bi?/i?+ (f (|) + 2||)B<0. 

It turns out that these dispersion relations (or at least their asymptotic form) 
are already known and have been used for the Navier-Stokes equation by J.-Y. 
Chemin in 0], and by D. Serre in |1J, 30 1 and B. Perthame in [Q. One of the 
interests of these estimates is that one may use them as an a priori estimate to 
control the behaviour for large times and build a global (in time) solution to the 
Cauchy problem. An equivalent remark (see J30)) is that it is possible to build 
a solution by a fixed-point method for a finite time (this is not in contradiction 
with T. Sideris' results |31 on non-existence, if the initial data is small in the 
correct sense) and that one may choose the rescaling ti-^ R{t) such that (for the 
complete rescaling as defined in Section 3 of course) the evolution with respect to 
the rescaled time holds only on a finite time interval 0<t<Too= L A~^ (t) dt. 
However, we are here interested only in the dispersion relations which were easily 
obtained by the mean of the second method of Section 4. These dispersion 
relations can be summarized as follows: 

Proposition 7.1 // (p,u) is a global classical solution of (IE) with 7> 1, then 
it satisfies the following dispersion relation 



d f n t 




R 


— R^ P 


n- 


X 


dt \ Jji^d 




R 



dx+^ f p\xfdx+-R'^f p''dx)<0 
R" Jw^ 7 Jtr'' J 



with (7 = min(2,(7-l)d) and t^R{t) such that R = R-^i+^'>, R{0) = Ra>0, 
R{0) = 0. 

8 Wigner and Schrodinger equations 

The relation between the Schrodinger equation, the Wigner equation and the 
Vlasov equation is now quite well understood. It has been the subject of a 
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considerable number of papers in the recent years: we mention ||l3, n3 as some 
of the most recent ones, and also 19, 23[ for the limit of the Schrodinger-Poisson 
to the Vlasov-Poisson system. Historically, the dispersion relations have been 
studied for the Schrodinger equation first and then adapted to the corresponding 
kinetic equation [^ ^, |l^ . The analysis of the dispersion relations in the kinetic 
framework came only after, but now seems to provide powerful tools to build 
new dispersion identities, cf. Q. 
Consider the Schrodinger equation 

If w is the Wigner transform of -ip, 

r % % 

w{t,x,v)= / e-''"y Tl>{t,x+ -y)->p{t,x - -y) dy , 
J^d z z 

it has to satisfy the Wigner equation 

1 
dtw + v-dx'w~—Q(V)w = 
n 

where the pseudo-differential operator Q{V) is defined by 
1 



eiv)f{x,v) 



i^T^f Jm" 



V{x+iy)-Vix-^y) 



TR'' 



In the semi-classical limit ?i-^0+, the operator Q{V) is formally expected to 
converge to —dxV- dy, and it is the purpose of many papers to justify this limit, 
cf. ll|,|i|,|l|,]2|. 

In this section we will only derive some dispersion identities according to the 
technique developed at the end of Section 4 and give some easy consequences 
of these estimates. 

We shall consider three cases: 

The linear case (L): y is a given fixed nonnegative potential which docs not 
depend on t and decays as |a;| ^ +oo. We will not go further into this case 
since the dispersion properties would depend on the local properties of V 
and x-dxV, but the co mpu tations are essentially the same as for the other 



cases up to Equation (8.1) 



The Poisson case (P): In case of the Schrodinger equation V obeys 

and in case of the Wigner equation V obeys 

—AV= / w{t,x,v)dv; 
Jtr'' 
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we consider only the electrostatic case. We shall state a result on the 
Wigner and the Schrodinger formulations of the problem, which clearly 
proves that this case can be handled in full generality with our methods. 
The estimates are slightly improved in dimension d — 3 and can obviously 
be generalized to any dimension rf>4. The results are new for d—2. 

The nonlinear case (NL): V is given by a power law 

In that case we consider only the nonlinear Schrodinger equation (NLS) 

(in the following, we shall only study the dcfocusing case e=— 1). This 
case is mentioned here to make the link with the pseudo-conformal meth- 
ods and to recover the pseudo-conformal law, which has been studied 
extensively. 

8.1 Wigner equation 

For the Wigner equation, we introduce as for the Vlasov-Poisson system the 
new velocity variable r]{t,x,v) = v— ^x and exactly as for the Vlasov-Poisson 
system, F{t,x,ri) — w{t,x,v) solves the rescaled Wigner equation (R'W): 

dtF + 7^-d,F-^x-d„F~^e{V)F+^(d,{xF)^dr,{vF)]^0. 

Again as for (R'VP), we compute the energy 

E{t)^ff F(\rj\^ + ^\x\^+av)dxdr] 

if (i> 3, and 

E{t)^^\ogR+ ff F(\r^\' + ^\x\' + av)dxdrj 

if (i = 2 in case (P). Here a is a coefficient which takes different values according 
to the case we consider: a — 2, 1, and -^ in case (L), (P) and (NL) respectively. 
The same computation as before provides 



dE 
~dt 



F{t,x,vi) 



^R , ,9 I d , R. ^R R \ , ,9 



dxdrj 



-a— j j V{t,x)do,{xF{t,x,r])) dxdri 
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(for (i>3 in case (P) — the case (P), d — 2 is similar up to the integrations by 
parts that are to be done with care) and we may define L{t) — B{t)E{t) and, as 
for the Vlasov-Poisson system, 



dL 






V{t,x)[ BF 



2^B 
R 



Bdx[xF) ) dxdrj 

R J 

F{t,x,ri)\ri\'^ dxdrj 



IR'^xIR'' 



«l 



d , R . R R , ^ 
dt^R^ + ^RR^'' 



M'^xTR.'' 



F{t,x,r])\x\ dxdrj. 



In the case of the coupling with the Poisson equation {d>2) , the conditions 
on L that are sufficient for it to be nonincreasing are exactly the same as for 
(R'VP) in the Poisson case (P): see Section 4. The detailed justifications of the 
computations for initial data 



mt=o,-)f 



w{t = 0,-,v)dveL\Wi'^) 



w 



are not given here, and we shall refer to [gj for a proof if (i> 3 in the context of 
the Schrodinger-Poisson system. 

Theorem 8.1 Assume that w is a solution of (WP) with M =\\w{t,-,-)\\ and 
thatt^R{t) is the solution of Equation (p^; i? + ei?i-'^ = 0, R{0) = 1, i?(0) = 
0. The function ti-^L{t) defined above for d>2 (with B = R'^^^) is decreasing 
for d = 2,2i, constant for d — A, and satisfies for any d> 2 



^^{d-A)RR'-^ 



R 



wdvdx . 



IJR.'^xM'^ 

In the plasma physics case e = —l, L is bounded from below and for d — 2,S, 



+ 00 



Ris)R''~%s) 



M'' X TR"^ 



R{s) 

V -—X 

Ris) 



w{s,x,v)dvdx I ds < +oo . 



However, the results on the dispersion for the Vlasov-Poisson system cannot be 
transposed straightforwardly because of the lack of positivity of w and one has 
to be very careful to recover the estimates given in p7[ for d — 3. In dimension 
d = 2, the situation is even worse because the boundedness of L from below is 
not obvious at all. In that sense, the Schrodinger formulation of the problem is 
more suitable. 

8.2 Schrodinger equation 

The Lyapunov function for the Schrodinger equation is easily found by simply 
considering the Wigner transform. However, it is interesting to realize how the 
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method of Section 4 applies directly. According to the Weyl quantification and 
the Wigner transform, the operator iTidx corresponds to the variable v. the 
change of variables rj = v— -^x therefore means that instead of iTidx we consider 
the new operator ihdx — ■4x: 

J. ^-ta ^ \^ -i-S^-ta f i-a^ 
(f)^^[iWx-—x)(j) = e « 2fi ihdxle li 2b 

For that purpose, we may consider the new wave function (j}{t,x) such that 

il;{t,x) = e'Tf^f^ (j)(t,x) . 
(j){t,x) solves the rescaled Schrodinger equation (R'S) 

^ndtcj,^-^h'A^+{V + ^\xf)^-'-^id^ + 2x-dx<t>). 

If we define the potential energy term by VF[0] = 21/|(/)p, M^[(/)] = y|0p, or 
M^[0] = -2-j-|(^|p+i in case (L), (P), or (NL) respectively, the corresponding en- 
ergy is given by 

E{t) = J^^(^n'\Wcf>\^ + ]¥[<!>] + ^\x\^\ci,\'yx 

if (i> 3, and 

i?W = j^^ (?i'|V0p + W^[0] + ||xn0|2^ dx + logi?(i)^(j||^J0pdx 

if d= 2 in case (P). We may then build the Lyapunov functional in the same way 
as for the solution of the Wigner equation. Going back to the original variables, 
we have to replace |V(?!)p by KV — i-^x)tp\'^. 

The Schrodinger-Poisson system and its asymptotics has been studied in 
[pi 0, 17 . More recently, a theory for L^ solutions corresponding to mixed 
quantum states has been established by F. Castella (see [§[§). In the case of a 
pure quantum state, J. L. Lopez and J. Soler in [pl| , 02] also gave detailed results 
on the asymptotic behaviour using a linear scaling approach in the continuation 
of the method developed S. Kamin and J. L. Vazquez. The main interest of our 
approach is that it gives a refined estimate for d — 3 and is adapted to the limit 
case d = 2 as well. 

Concerning the notion of solution we may assume that it is as smooth as 
desired and refer to p, H] for minimal requirements (estimates for weak solutions 
are built using approximating smooth solutions). 
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Theorem 8.2 Assume that d>2 and consider a solution of the Schrodinger- 
Poisson system. With the above notation 



L{t) = R'^-\t) 
for d== 3,4, and 

m= 






dx + R^-'it) / V\'il;\'dx + 



1 



RHt) 



\x\^\^\^dx 






-v\M' 



'-ri^p)d.+M%miiogi?(,) 



R'it) 



2n 



for d—2, is decreasing for d—2, 3 and constant for d — A if t^-^ R{t) is a so- 
lution of R = R^^'^, R{0) — 1, i?(0) = 0. As a consequence n{t,x) = \'ip{t,x)\'^ is 
decreasing: there exists a constant C > such that 



ll-(i,-)IL.(H^)<c.7?'^(§-i).i?'^(i-^)(i- 



(8.1) 



for any pe [2,^] if d^3, 4 and liminff^+oo ||"(i,-)llLp(]R2) =0 for pe]2,+oo[ 
ifd = 2. 

Note that for d — 3, p—10/3, we recover the same exponents as for the 
Vlasov-Poisson system. For d = 2, exactly the same estmiate as in the proof of 



Theorem 4.2 holds 



liminf 

t^ + oo 



]R2 



(V~z— .)^ 



dx^Q. 



The crucial ingredient in the proof of this theorem is the following interpolation 
lemma (se e H , Q] and ||l^, Cor. 5.5]) which plays a role similar to the one of 
Equation ( ^.l| ) for the Vlasov-Poisson system: 

Lemma 8.3 Assume that d>i. There exists a constant C>0 depending only 
on d such that, for any u£H^{iR ) such that xi-^xu{x) belongs to L^{]R ), 

M\L^^n'^~,<C\\u\\l,^^,^\\ix + itV)u\\l-^^,yt-^'--^ 
foranype[2,^^],a^l{l-^). 



The proof of Lemma B.3 is easily established using the Gidas-Nirenberg inequal- 
ity 

hllLP(R^) < [C(d)]'""ll«llL.(R.)l|V^||l-JL 

where C{d) is the Sobolev constant corresponding to the injection of H^{Si ) 
into L"^^^ (JR ), and the decomposition u — pe^"^ which holds at least for smooth 
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enough functions (the conclusion holds by a density argument). We may then 
write 

\\{x + itV)u\\l2=t'^ / \Vpfdx+ / \xp + tpVipf dx>t'^\\V\u\\\l2 , 

which proves the interpolation results. 



Proof of Theorem |8.2| . One has to replace 1/t by R/R: for d~3 or 4, 
we may refer to |0 for the proof of Equation (|^), where it is done in the 
case R — t. For d — 2, the argument is similar, the main step being the proof of 
the boundedness of L which goes exactly as in the Vlasov-Poisson case (see the 



proof of Theorem ^) . D 



We conclude this section by considering the case of the Nonlinear Schrodinger 
equation which allows us to make an explicit link with the pseudo-conformal law. 
If W[(j)] = -rilf/'l^^"^, a direct computation gives 

dE d{p-l)R f^^^, ^^ , „^^2/",^,,2, f d ,R, „^-R\ /", ,2, ,,2 , 

^ ^ 'W[(j)\dx-2 — ?i\V(t)fdx+[—{ — ) + 2 — —]\xf\(t)fdx 



dt 2R J ' ' R J \dt'R/ RR 

and L{t)^B{t)E{t) is decreasing if B{t)^Ri{t), g = min((p- l)d/2,2), R = 
1/i?^"''^, i?(0) = l, i?(0) = 0. In the next result we are again not interested in 
the weakest notion of solution and assume that the solution is global in t and 
as smooth and sufficiently decreasing at spatial infinity as necessary to justify 
any integration by parts in the computations. 

Theorem 8.4 Assume that d>2 and consider a global solution of the Nonlinear 
Schrodinger equation (NLS) 

Then with the above notation 

Lit)^R^it)|(^\iV-^^x)i;\^ + -^\yJf+'yx+J^J\x\'\^\'dx (8.2) 
is decreasing. 



Decay estimates can of course be deduced from Lemma B.5 as for the 
Schrodinger-Poisson system. The details of the computations for the proof of 
Theorem ^.41 are left to the reader. 

A simple method to understand the pseudo-conformal law is simply to look 
for a pseudo-conformal invariance of the equation, i.e. a transformation which 
leaves the equation invariant. If u{t,x) is a solution of (NLS) in the focusing or 
in the defocusing case (e = — 1): 

1 2 
ihdtu — — h Au + e\u\P u, 
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let us look for a function ti-^ {R{t),T{t),uj{t)) such that (r,^) t— >w(t,^) given by 



is a solution of (NLS) too for some a G ]R. It turns out that this is possible only 
in the case p— 1 = | (critical case), and in that case, ti-> (i?(i),r(t),cj(i)) solves 
the system 

dr 1 dR „ „ duj „ o 

dt R^' dt ' rfi 

The solution is given by 

^W- , /? , , R{t) = Roil + 2iJot), r(t)^ [ +ro. 



This transformation can be found in |gj| (see also |18 for instance). The con- 
servation of the energy after rescaling (conservation of the energy for v) gives 
the following conservation law for u: 

— (R'^{t)[ \\/u{t,x)^iuj{t)xu{t,x)\^dx--^ [ \u{t,x)\^^''+^Ux^=0. 
dt\ jRd d + 2 7]Rd J 

(8.3) 

This expression clearly corresponds to the case q = 2 — {p—l)d/2, and the 

pseudo-conformal law is nothing else than the expression of dL/dt where L 

is given by Equation ( ^.2[ ). As we already noticed already several times, one 

may replace uj(t) and R(t) by their equivalents as t—^+oo, which is the same 

as considering the singular solution corresponding to the limit luq -^ -f-oo and 

i?oWo -^ Ij and recover instead of Equation ( |8.3| ) the more classical form for the 

conformal invariance law: 



d 
di 



w 



V u{t,x) ~ i—u{t,x) 



de 
dx — 



d + 2 



f u{t,x)\i'-''+^Ux)^0. 
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